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and University of California, Berkeley 

Given a two-level regular fractional factorial design of resolu- 
tion IV, the method of doubling produces another design of resolu- 
tion IV which doubles both the run size and the number of factors 
of the initial design. On the other hand, the projection of a design of 
resolution IV onto a subset of factors is of resolution IV or higher. 
Recent work in the literature of projective geometry essentially de- 
termines the structures of all regular designs of resolution IV with 
n > N/4 + 1 in terms of doubling and projection, where N is the run 
size and n is the number of factors. These results imply that, for 
instance, all regular designs of resolution IV with 5N/1Q < n < N/2 
must be projections of the regular design of resolution IV with N/2 
factors. We show that, for 9A r /32 < n < 5N/16, all minimum aberra- 
tion designs are projections of the design with 5iV/16 factors which 
is constructed by repeatedly doubling the 2 5_1 design defined by 
I — ABCDE. To prove this result, we also derive some properties of 
doubling, including an identity that relates the wordlength pattern 
of a design to that of its double and a result that does the same for 
the alias patterns of two-factor interactions. 



1. Introduction. Doubling is a simple but powerful method of construct- 
ing two-level fractional factorial designs, in particular, those of resolution IV. 
Suppose X is an N X n matrix with two distinct entries, 1 and —1. Then 
the double of X, denoted by D(X), is the 2N x 2n matrix 
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that is 



D(X) 



1 1 
1 -1 



where (g) is the Kronecker product. Suppose X defines an A-run design for 
n two- level factors, where the two levels are denoted by 1 and —1, each 
column of X corresponds to a factor and each row of X defines a factor- 
level combination. Then -D(X) defines a design which doubles both the run 
size and the number of factors of X. Such a method was used by Plackett 
and Burman [10] in their classical paper on orthogonal main-effect plans. 

An N xm submatrix of X, where m<n, is called a projection of X (onto 
m factors). Equivalently, such a design can be obtained by deleting n — m 
columns (factors) from X. We allow the possibility m = n, so that a design 
is considered as its own projection onto all columns. If X is of resolution IV, 
then all projections of X are of resolution IV or higher. A less obvious fact 
is that if X is of resolution IV, then -D(X) is also of resolution IV. Thus, 
starting with a design of resolution IV, the combined operation of doubling 
followed by projection yields a design of resolution IV or higher. 

Some recent results in the literature of finite projective geometry [2, 3, 8] 
essentially characterize, in terms of doubling and projection, the structures 
of regular designs of resolution IV with n > N/4 + 1. These elegant and 
difficult results have important implications in statistical design, but are 
not easily accessible to statisticians since they were stated in the language 
of projective geometry. One purpose of this paper is to review some of these 
results and rephrase them in design language. This should be useful to people 
interested in the statistical design of experiments. In the past few years we 
have encountered several occasions where some special cases of these results 
were either re-discovered by statisticians or would have been very helpful if 
they had been become more widely known in the statistical literature. 

We also present some new results, including an identity that relates the 
wordlength pattern of a design to that of its double and a result that does 
the same for the alias patterns of two- factor interactions. These results and 
some other basic properties of doubling are presented in Section 2. Char- 
acterization of regular designs of resolution IV with n > A/4 + 1 in terms 
of doubling and projection is discussed in Section 3. Section 4 shows how 
the results in Sections 2 and 3 can be used to investigate certain minimum 
aberration designs. In particular, we show that, for 9A/32 < n < 5 A/16, a 
minimum aberration design can be obtained by deleting factors from the 
minimum aberration design with 5N/16 factors which can be constructed 
by repeatedly doubling the 2 5_1 design defined by / = ABCDE. 

We conclude this section with a review of some basic terminology and 
notation. A factorial design is called an orthogonal array of strength t if, in 
all projections onto t factors, all factor-level combinations appear the same 
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number of times. Regular fractional factorial designs are those which can be 
constructed by using defining relations and are discussed in many textbooks 
on experimental design; see, for example, [11]. Each interaction that appears 
in the defining relation is called a defining word, and the resolution of a 
regular design is defined as the length of the shortest defining word. It is 
well known that a regular design of resolution R is an orthogonal array of 
strength R — l. For each positive integer i, let B{ be the number of defining 
words of length i. Then the resolution is equal to the smallest i such that 
Bi > 0. The sequence {B\ , B2, ■ ■ ■ , B n ) is called the wordlength pattern of the 
design. The minimum aberration criterion introduced by Fries and Hunter 
[9] chooses a design by sequentially minimizing B\, B2, B3, — 

Throughout this paper N and n denote the run size and the number of 
factors, respectively. We shall restrict consideration to two-level designs only. 
Under a regular two-level design, N must be a power of 2, say, N = 2 n ~ p . 
This design is a ^-fraction of a complete 2 n factorial, and is usually re- 
ferred to as a 2 n ~ p design. It is well known that a regular two-level design 
of resolution III must have n < N — 1 , and a regular two- level design of res- 
olution IV must have n < N/2. A design of resolution III is called saturated 
if n = N — 1, and a design of resolution IV is called saturated if n = N/2. 

2. Some basic properties of doubling. We first note that if X is a 
Hadamard matrix of order N, then D(X.) is a Hadamard matrix of order 
2N. In fact, this was what Plackett and Burman used in their 1946 classic 
paper. The following properties can easily be established: 

Theorem 2.1. If X is an orthogonal array of strength two, then -D(X) 
is also an orthogonal array of strength two. Likewise, if X is an orthogonal 
array of strength three, then -D(X) is an orthogonal array of strength three. 

Theorem 2.1 has no counterpart for designs of higher strength. In fact, 
for two columns a and b of X, -D(X) must have four columns of the form 

lab-a-bl' whose componentwise product has all the entries equal to 1. 
Therefore, D(X.) cannot have strength higher than three. For regular de- 
signs, these four columns would lead to a defining word of length four. 

Another important elementary fact is that saturated regular designs of 
resolution III are unique (up to isomorphism). Such a design of size 2 k can 
be obtained by deleting the first column of 
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In other words, saturated regular designs of resolution III can be obtained 
by deleting a column of l's after successively doubling the 2x2 matrix [J ]]. 
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Since all regular designs of resolution III or higher can be constructed by 
deleting a subset of columns from a saturated regular design of resolution III 
(or, equivalently, by projecting the saturated regular design of resolution III 
onto the complementary set of columns), all two- level regular designs of 
resolution III or higher can be constructed by the operation of doubling fol- 
lowed by deletion (or projection). These observations also imply that if X is 
a regular design, then -D(X) is regular. In other words, in Theorem 2.1, if X 
is a regular design of resolution III (resp. IV), then -D(X) is also a regular 
design of resolution III (resp. IV). However, by the comment following The- 
orem 2.1, if X is a regular design of resolution higher than IV, then -D(X) 
is a regular design of resolution IV only. 

Saturated regular designs of resolution IV are also unique (up to iso- 
morphism). One way of constructing such designs is to foldover saturated 
regular designs of resolution III. We would like to present another neat and 
compact method of construction using doubling. Consider the 2 2 complete 
factorial 



whose number of factors is precisely half of the run size. Successively dou- 
bling the 2 2 complete factorial results in designs whose number of factors 
is half of the run size, and by the discussion at the end of the previous 
paragraph, are regular designs of resolution IV. This implies that saturated 
regular designs of resolution IV can be obtained by successively doubling 
the 2 2 complete factorial. 

As mentioned earlier, the designs obtained by doubling those of resolu- 
tion III or higher are of resolution III or IV. To study statistical properties of 
such designs, it is important to consider the alias pattern of two-factor inter- 
actions. Suppose N = 2 n ~ p . Among the 2" — 1 factorial effects, 2 P — 1 appear 
in the defining relation. The rest are divided into g = 2 n ~ p — 1 alias sets, each 
of size 2 P . Without loss of generality, assume that the first / = 2 n ~ p — 1 — n 
of these alias sets does not contain main effects. For each 1< i < g, let rrii be 
the number of two- factor interactions in the ith alias set. Cheng, Steinberg 
and Sun [7] showed that 
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(2.1) 




and 



(2.2) 




DOUBLING AND PROJECTION 



5 



Furthermore, estimation capacity, a measure of model robustness discussed 
in [7], can be expressed in terms of the m^s. The following result relates the 
rrii values of a design to those of its double. 

Theorem 2.2. Suppose a regular design X has g alias sets, the first 
f of which do not contain main effects, and for 1 <i < g, rrii is the number 
of two- factor interactions in the ith alias set. Let the corresponding num- 
bers of D(X.) be g*, f* and m*, respectively. Then g* = 2g + 1, f* =2f + 1, 
ml = m\ = 2mi, mg = m\ = 2rri2, ■ ■ ■ , m* 2 ^_ x = = 2m/, m<2j +1 = n, and 
m 2/+2 = m 2/+3 = 2m /+i; ■■■,m**_ 1 = m* g * = 2m g . 

Proof. Suppose X is a 2 n ~ p design. Then D(X) is a 2 2 ™-(™+p- 1 ) design, 
and g = 2 n ~P - 1, / = 2 n ~P -l-n,g* = 2 2n ~( n +P~ 1 ) - 1, /* = 2 2 ™-(™+p- 1 ) - 
2n - 1. The relations /* = 2/ + 1 and g* = 2g + 1 follow. 

To each factor in X, say A, there correspond two factors in D(X). Suppose 
the column of X corresponding to A is a. We shall denote the factor in 
-D(X) corresponding to the column [ a ] by A + and the factor corresponding 

to [_ a ] by A~ . Then it is easy to see that, for any two factors A and B 
in X, under -D(X) A + B + and A~B~ are aliased ,A + B~ and A~ B + are 
aliased, and A + A~ and B + B~~ are aliased. Consequently, if the two-factor 
interactions AB and CD are aliased under X, then A + B + , A~B~, C + D + 
and C~D- are aliased under D(X), and A+B~, A-B+,C + D~ and C~D+ 
are also aliased. If the main effect A and two- factor interaction CD are 
aliased under X, then A + , C + D + and C~D~ are aliased under L>(X), and 
A~ , C + D~ and C~D + are also aliased. Using these facts, it can be seen 
that each alias set of two-factor interactions under X determines two alias 
sets of two- factor interactions under D(X), each of which is twice as large as 
the original alias set under X. Furthermore, all the n two-factor interactions 
of the form A + A~ constitute another alias set which does not contain main 
effects. □ 

We have the following relationship between the wordlength pattern of X 
and that of D(K). 
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Theorem 2.3. Let Bk and be the number of defining words of length 
k o/X and Z)(X), respectively. Then 

( min[(fc-l)/2,n] , . 



=o 



min[fc/2-l,n] . / \ 

E 5 fe _ 2s .( n -^ +2s )2^- i + ^j, 

i//c is ei>en and is even; 

min[fc/2-l,n] / n _ /j _|_ 2s \ 

E -Sfc-2s " ( ) 2 fc ~ 2s_1 , if k is even and k/2 is odd. 



s=0 



Proof. If Ai t ■ ■ ■ Ai k is a defining word of X, then A?* ■ ■ ■ A^ k , where each 
ji is a + or — , is a defining word of D(X.) as long as an even number of the ji 's 
are — 's. When s is even, for any s distinct factors At, , .... A , A? A~ ■ ■ ■ Af Ay 
is also a defining word of Z)(X). In general, a defining word of -D(X) is of the 
form Af Ay ■ ■ ■ At Ay A j y +1 ■ ■ ■ A J - 3+t , where As, , . . . , At , . are distinct factors 
of X, Ai s+1 ■ ■ ■ Ai g+t is a defining word of X, and an even (resp. odd) number 
of the ?Vs are — 's when s is even (resp. odd). If At Ay ■ ■ ■ At Ay A J - 3+1 • • ■ A 3 y +t 

> r ' II '1 l s Is «s + l *s + i 

is of length k, then t = k — 2s. This implies that < s < min[/c/2,n]. On 
the other hand, when s is odd, we must have t > 1. Therefore, we have 
s < min[(/c — l)/2,n] when is odd, s < min[fc/2,n] when /c is even and /c/2 
is even, and s < min[fc/2 — l,n] when k is even and k/2 is odd. The theorem 
then follows from the fact that, for given s, there are 2* _1 = 2 fc_2s_1 ways 
to choose (js+i, ■ ■ ■ ,j s +t) if * > 1, B k ^ 2s ways to choose A is+1 , . . . , A is+t and 
( n ~^+ 2s ) ways to choose A^, . . . ,Ai 3 . Note that the last term in the case 
where k/2 is even corresponds to s = k/2. In this case, there are (fe)ways to 
choose At Ay ■■ -At Ay . □ 

In particular, for designs of resolution III or higher we have B% = AB^, and 
B\ = 8B4 + (2). These identities can also be derived by using (2.1), (2.2) and 
Theorem 2.2. 

The following is an immediate consequence of Theorem 2.3. 

Corollary 2.4. Given two regular designs Xi and X2, Xi has less 
aberration than X2 if and only if D(X.\) has less aberration than D(X.2). 

3. Maximal designs of resolution IV. One major difference between de- 
signs of resolution III and IV is that all regular designs of resolution III can 
be constructed by deleting factors from saturated regular designs of resolu- 
tion III. On the other hand, not all regular designs of resolution IV can be 
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obtained by deleting factors from saturated regular designs of resolution IV: 
while all designs obtained by deleting factors from saturated regular designs 
of resolution IV are the so-called even designs which have no defining words 
of odd lengths, there are designs of resolution IV which are not even designs. 

We say that a regular design of resolution IV or higher is maximal if 
its resolution reduces to three whenever an extra factor is added. Clearly, 
the saturated regular designs of resolution IV are maximal. If a design is 
not maximal, then at least one factor can be added so that the design is 
still of resolution IV or higher. One can keep adding factors until it becomes 
maximal. Therefore, if a regular design of resolution IV is not maximal, then 
it can be obtained by deleting factors from a maximal design. Because of 
the importance of this fact, we state it formally as a proposition: 

Proposition 3.1. Every regular design of resolution IV is a projection 
of a certain maximal regular design of resolution IV or higher. 

We can also define maximal regular designs of resolution III. It turns 
out that, for a given run size, there is only one maximal regular design of 
resolution III: the saturated regular design of resolution III. As mentioned 
earlier, saturated regular designs of resolution IV are maximal, but they are 
not the only maximal designs of resolution IV. 

Theoretically speaking, for a fixed number of runs, if we can determine 
all the maximal designs of resolution IV or higher, then all regular designs 
of resolution IV can be constructed by projecting the maximal designs onto 
subsets of factors. Recently some significant progress in the determination of 
such maximal designs has been made in the literature of finite projective ge- 
ometry. Note that maximal designs of resolution IV or higher are equivalent 
to maximal caps in a finite projective geometry. 

We first state a simple but interesting characterization of maximal regular 
designs of resolution IV. 

Theorem 3.2. A regular design of resolution IV is maximal if and only 
if mi > for all i = 1, . . . , /. 

This result was first stated in the coding-theoretic language; see, for ex- 
ample, [2]. Chen and Cheng [5] rephrased it in the above form. They also 
defined the notion of estimation index. Another way to state the result in 
Theorem 3.2 is that a regular design of resolution IV is maximal if and only 
if its estimation index is equal to 2. Since we can estimate one effect from 
each alias set assuming that the other effects in the same alias set are negli- 
gible, Theorem 3.2 says that a regular design of resolution IV is maximal if 
and only if all the available degrees of freedom can be used to estimate main 
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effects and two-factor interactions. Such designs are said to be second-order 
saturated by Block and Mee [1]. 

The following result, whose geometric version can also be found in [2], 
reveals the crucial role played by the method of doubling in constructing 
designs of resolution IV. 

Theorem 3.3. Let X be a regular design of resolution IV or higher. 
Then X is maximal if and only if D(X.) is maximal. 

One can see that Theorem 3.3 follows immediately from Theorem 2.2 and 
Theorem 3.3. 

The following two key results essentially determine the structures of reg- 
ular resolution IV designs with N/4 + 1 < n < N/2. 

Theorem 3.4 ([2, 8]). Every maximal regular design of resolution IV 
with N/4 + 2 <n< N/2 can be obtained by doubling a maximal regular design 
of resolution IV or higher. 

Theorem 3.5 ([3]). For each N = 2 k with k>4 there exists at least one 
maximal regular design of resolution IV or higher with n = N/4 + 1. 

One family of maximal designs with n = N/4 + 1 can be found in Tang, 
Ma, Ingram and Wang's [13] study of designs with maximum number of 
clear two-factor interactions. 

A complete search shows that there are two maximal regular 16-run de- 
signs of resolution IV or higher. One is the saturated 2 8-4 design of reso- 
lution IV, and the other is the 2 5_1 design defined by I = ABCDE, whose 
existence is assured by Theorem 3.5. Repeatedly doubling the former yields 
all larger saturated regular designs of resolution IV, while successively dou- 
bling the latter (a design of resolution V) leads to a family of maximal 
designs of resolution IV with n = 5JV/16, for N = 2 k , k > 5. Since there 
are no 16-run maximal regular designs of resolution IV with 5 < n < 8, by 
Theorem 3.4, for all N = 2 k with k > 5, there are no maximal designs of 
resolution IV with 5iV/16 < n < N/2. This leads to the following important 
conclusion. 

COROLLARY 3.6. For 5N/16 < n < N/2, a regular design of resolu- 
tion IV must be obtained by deleting columns from saturated designs of res- 
olution IV. All such designs are even designs. 

Doubling the two maximal 16-run regular designs of resolution IV, we 
obtain the saturated 2 16-11 design of resolution IV and a 2 10 ~ 5 design, both 
of which are maximal. By Theorem 3.5 there exists at least one maximal 
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regular 32-run design of resolution IV with n = 9. In fact, there is exactly 
one such design. Therefore, for 9 < n < 16 there are exactly three 32-run 
maximal designs of resolution IV: a 2 16-11 , a 2 10-5 and a 2 9-4 . The last 
one and its repeated doubles constitute a family of maximal regular designs 
of resolution IV with n = 97V/32, for N = 2 k , k > 5. A gain, there are no 
maximal regular designs of resolution IV with 9iV/32 < n < 5iV/16. Thus, 
for the n's in this range, a regular design of resolution IV must be obtained 
by deleting columns from either the saturated design of resolution IV or the 
maximal regular design of resolution IV with n = 5N/16. 

Doubling the three maximal 32-run designs of resolution IV, we obtain 
the saturated 2 32-26 design of resolution IV, a 2 20 ~ 14 design and a 2 18 ~ 12 
design, all of which are maximal. By Theorem 3.5 there exists at least one 
maximal regular 64-run design of resolution IV with n = 17. Block and Mee's 
[1] complete search shows that there are five such designs. Therefore, for 
17 < n < 32 there are eight 64-run maximal regular designs of resolution IV: 
a 2 32 " 26 , a 2 20 " 14 , a 2 18 ~ 12 and five 2 17 " 11 . One can also conclude that there 
are exactly five maximal regular designs of resolution IV with n = 17iV764, 
for all N = 2 k with k > 6. 

Now it is clear that if N = 2 k , k > 4, then for n > N/4 + 1 a maximal 
regular design of resolution IV or higher must have 

(3.1) n G {N/2, 5N/W, 9N/32, 17A/64, 33A/128, . . . }. 

Conversely, for each integer n = (2* + l)N/2 l+2 , there exists at least one 
maximal regular iV-run design of resolution IV or higher with n factors. A 
maximal regular design of resolution IV or higher with n = (2* + \)N/2 l+2 
and N = 2 k , k> i + 2, can be constructed by repeatedly doubling a maximal 
regular 2 J+2 -run design of resolution IV or higher with 2 % + 1 factors. 

4. Some results on minimum aberration designs. By the discussion in 
the previous section, there are maximal regular designs of resolution IV 
with n = N/2, 5N/16, 9A/32, 17A764, 33A7128, .... Those with n = N/2, 
the saturated regular designs of resolution IV, are known to have minimum 
aberration. For 5iV/16 < n < N/2, minimum aberration designs (in fact, all 
regular designs of resolution IV) must be projections of the saturated regular 
design of resolution IV. Butler [4] addressed the issue of deleting factors from 
the saturated regular design of resolution IV so that the resulting design 
has minimum aberration. This is reminiscent of the complementary design 
theory of Chen and Hedayat [6], Tang and Wu [14] and Suen, Chen and Wu 
[12] that deals with how to find a set of factors so that its complement in 
the saturated regular design of resolution III has minimum aberration. 

In an unpublished work, N. A. Butler found the minimum aberration 
design with n = 5N/16 which, using the terminology in this paper, is the 
maximal regular design of resolution IV with n = 5-/V/16. In this section 
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we shall show that, for 9iV/32 <n< 5N/16, the minimum aberration de- 
signs are projections of the maximal regular design of resolution IV with 
n = 5N/16. Thus, although for n > 9iV/32 minimum aberration designs are 
projections of the maximal regular design of resolution IV with either N/2 
or 5iV/16 factors, the first two in (3.1), the pattern breaks down at 9iV/32. 
Even the maximal design of resolution IV with n = 9iV/32 itself does not 
have minimum aberration. 

Before proceeding to the proof of this result, we shall explore a bit more 
the alias pattern of two-factor interactions under the maximal design of res- 
olution IV with n = 5N/16. Suppose N = 16 • 2 , where t > 1, and let X* be 
the maximal design of resolution IV with 5 • 2* factors. Then X* can be ob- 
tained by doubling the 2 5_1 design defined by / = ABODE t times. Suppose 
a, b, c, d and e are the five columns of this 2 5_1 design corresponding to fac- 
tors A,B,C,D and E, respectively. Then each of a, b,c,d and e generates 
2* columns of X*. For example, each of the 2* columns of X* generated by 
a takes the form xt <8) • • • <8) xi Cg> a, where Xj = [1, 1] T or [1, — l] T . We shall 
denote the corresponding factor of X* by A 3 , where j = (ji, ■ ■ ■ ,jt), with 
ji = 1 if Xi = [1, 1] T and ji = -1 if x, = [1, -1] T . Notation such as B 3 , C j , D } 
and E } is similarly defined. The 5 ■ 2* factors of X* are thus partitioned into 
five groups each of size 2*. 

Any two of the five factors A,B,C,D and E, say X and Y, generate 
2* • 2* two-factor interactions of the form X*Y 3 , where i and j are 1 x t 
vectors with entries 1 or —1. These interactions, (2)2* • 2* = 10 • 2* • 2* in 
total, are called between- group two- factor interactions. Those of the form 
X*X } with i 7^ j, 5(2) = 5 • 2*~ 1 (2 i — 1) in total, are called within-group 
two-factor interactions. 

Since the 2 5_1 design defined by / = ABODE is of resolution V, there is 
exactly one two-factor interaction in each of its ten alias sets not containing 
main effects. By applying Theorem 2.2 t times, we see that, under X*, there 
are 10 ■ 2* + (2* — 1) alias sets not containing main effects, 10 • 2 t of which each 
containing 2* two-factor interactions, and 2* — 1 of which each containing 
5 • 2* _1 two- factor interactions. 

It can be seen that the 2 t ■ 2 t between-group two-factor interactions arising 
from the same pair (X, Y) are distributed evenly in 2* alias sets of size 2 t . 
Each of these 2* alias sets is of the form {X'yJ : i j = k}, where k is a 1 x t 
vector with entries 1 or —1 and iQj is the componentwise product of i and j. 
We denote this alias set by XY^. The ten possible pairs of X and Y account 
for the 10 • 2* alias sets of size 2* mentioned in the previous paragraph. On 
the other hand, the 2* _1 (2* — 1) within-group two-factor interactions X 1 X i 
arising from the same X are distributed evenly in the remaining 2* — 1 alias 
sets. Each of these alias sets consists of the 5 • 2 t ~ 1 interactions A 1 A 1 , B X B S , 
C'CJ, D i D i , E l Ei with i©j = k. (Note that X i X i is the same as X'^X 1 and 
k^ 1, where 1 is the vector of l's.) We denote this alias set by W^. 
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One key property that is important for the proof is that each of the 2 
factors of the form X 1 appears in exactly one of the 2* two-factor interactions 
in each alias set XY^. As a consequence, if u factors generated by X are 
deleted from X*, then the number of two- factor interactions in each of these 
2 t alias sets is reduced by u. In this case, since each factor of the 2 5_1 design 
can be coupled with four other factors, the number of two- factor interactions 
in 4 • 2* of the 10 • 2* alias sets of size 2* is reduced to 2* — u; that in each of 
the other 6 ■ 2* alias sets remains to be 2* . 

It can also be seen that each of the 5 • 2* factors appears in exactly one 
within- group two- factor interaction in each W^. 

Now we are ready to prove the following theorem. 

Theorem 4.1. For any N = 2 k , k>5, and 9A732 <n< 5N/W, the 
minimum aberration design must be a projection of the design constructed 
by repeatedly doubling the 2 5 " 1 design defined by I = ABCDE. 

Proof. Let N = 16-2* and u = 5N/16 - n. Then since 5N/16 - 9iV/32 = 
2 t ~ 1 , we have 0<u<2* _1 . 

The minimum aberration design with N = 32 and n = 9 is known to have 
at least one zero among the m;'s, 1 < i < f ([7], page 91), and therefore is 
not maximal. In fact, it is obtained by deleting one factor from the maximal 
design with ten factors. By repeatedly applying Corollary 2.4, we see that, 
for each N = 2 k with k > 5, the maximal regular design of resolution IV 
with 9N/32 factors does not have minimum aberration. Now since there are 
only two maximal regular designs of resolution IV with more than 9iV/32 
factors, one with N/2 factors and the other with 5N/16 factors, it is enough 
to show that, for N = 16 • 2* with t>2, there is at least one (5-/V/16 — u)- 
factor projection of the maximal regular design of resolution IV with 5iV/16 
factors that has less aberration than all projections of the saturated regular 
design of resolution IV. As before, let X* be the maximal regular design of 
resolution IV with 5N/16 factors. 

Now consider a design X obtained by deleting from X*« factors that are 
generated by the same factor A of the 2 5_1 design defined by I = ABCDE. 
We shall show that X has fewer defining words of length four or, equivalently, 
by (2.2), a smaller value of X)f=i m i than all projections of the saturated 
regular design of resolution IV. 

By the observations preceding the statement of the current theorem, under 
X, 4 • 2 t of the 10 • 2* alias sets of between-group two-factor interactions have 
mi equal to 2 t — u, and the other 6 • 2* alias sets have rri{ equal to 2 t . Now 
we provide an upper bound on the sum of squares of the m, values over the 
alias sets of within-group two- factor interactions. 

Under X* , in each of the 2* — 1 alias sets of within-group two- factor inter- 
actions, 2* _1 of the 5 • 2* _1 two-factor interactions involve factors generated 
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by A. Thus, when only factors generated by A are deleted, for each of these 
alias sets, the resulting mi satisfies 4 • 2 4 " 1 < m; < 5 • 2* . Consequently, an 
upper bound on the sum of squares of the mi values over these alias sets can 
be obtained by making as many of the m; values equal to 5 • 2 t_1 or 4 • 2* _1 
as possible. 

When u factors generated by A are deleted, Q) + u(2 l — u) within-group 
two-factor interactions are also deleted. Write 

u(u + l)/2 = a-2 t ~ 1 + b, 

where a and b are nonnegative integers such that b < 2* . Then 

(4.1) W +u{2 1 -u) = {2u-a)2 t ~ l -b. 

By the observation two paragraphs above, an upper bound on the sum of 
squares of the rrii values over the 2* — 1 alias sets of within-group two-factor 
interactions can be obtained by assuming that 2u — a — \ of the rrii values are 
equal to 4 • 2 t_1 , one is equal to 4 • 2 t ~ 1 + b, and the remaining 2* — 2u + a — 1 
values are equal to 5 • 2*" 1 . Combining this with the rrii values for the alias 
sets of between-group interactions obtained earlier, we conclude that if u 
factors generated by A are deleted from X*, then 



/ 

^m 2 < 4 • 2* • (2* - uf + 6 • 2* • (2*) 2 + (4 • 2*" 1 + bf 
(4.2) 1=1 

+ (2* - 2u + a - 1)(5 ■ 2*" 1 ) 2 + (2« - a - 1)(4 • 2*- 1 ) 2 . 

On the other hand, a saturated regular design of resolution IV has all the 
( 5 ' 2 2 _n ) two- factor interactions in N/2 — 1 = 2* +3 — 1 alias sets. Thus, a design 
obtained by deleting factors from a saturated regular design of resolution IV 
can have at most 2 t+3 — 1 nonzero m/s. For such a design, 



/ 

(4.3) £m 2 > 
=1 



r;-T 2 

2*+ 3 - 1 



• (2 i+3 - 1) 



It is sufficient to show that the right side of (4.2) is less than that of (4.3). 
This can be verified by tedious calculations using (4.1) and the assumption 
that < b < 2 t_1 . The details are omitted. □ 

Note that the case u = provides an alternative proof of Butler's result 
on the optimality of maximal designs of resolution IV with n = 5N/16. 

Theorem 4.1 leaves open the issue of which projection of the maximal 
design of resolution IV with n = 5N/16 has minimum aberration. A com- 
plementary design theory in the same spirit as that for saturated regular 
designs of resolution III and IV needs to be developed. 



DOUBLING AND PROJECTION 
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Suppose u factors are deleted from X*. For X = A,B,C,D,E, let nx be 
the number of factors deleted from those generated by X. Among the nxny 
pairs (i,j) where X 1 and Y 3 are deleted, let n XY be the total number such 
that i j = k; this is the number of between-group two-factor interactions 
formed by the deleted X 1 and y j 's that belong to the alias set XY^. Simi- 
larly, among the ( n ^) ordered pairs (i,j) where X 1 and X 3 are deleted, let 
n x x be the total number such that i j = k; this is the number of within- 
group two- factor interactions formed by the deleted A"s that belong to the 
alias set W^. Then the number of two- factor interactions in XY^ is reduced 
by nx + ny — n|y, and the number of two-factor interactions in is re- 
duced by ua + riB + nc + n D + n E - n\ A - n\ B - n^j C - n DD - n\ E = 

u ~ n \ A~ n \ B~ n \; c~ n D D~ n E E- Thus, a minimum aberration projec- 
tion of X* must minimize 

X,Y k 

+ E( 5x - u + n %A + + nc,c + ^b.D + n E , E ) 2 , 

where, in the first term, the first sum is over the ten possible (X, Y)'s and 
the second sum is over all the 1 x t vectors of l's and — l's. While this can 
be solved without difficulty for small u's, more general results need to be 
developed. We expect an optimal strategy to delete the factors one at a 
time alternately from the five groups of factors generated by A,B,C,D,E, 
while making the two-factor interactions formed by the deleted factors as 
uniformly distributed among the alias sets as possible. 

Acknowledgments. We would like to thank the referees for their helpful 
comments. 
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